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K-essence and the coincidence problem 
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K-essence has been proposed as a possible means of explaining the coincidence problem of the 
Universe beginning to accelerate only at the present epoch. We carry out a comprehensive dynamical 
systems analysis of the k-essence models given so far in the literature. We numerically study the 
basin of attraction of the tracker solutions and we highlight the behaviour of the field close to sound 
speed divergences. We find that, when written in terms of parameters with a simple dynamical 
interpretation, the basins of attraction represent only a small region of the phase space. 
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I. INTRODUCTION 



The search for an explanation for the observed accel- 
eration of the Universe 0| is one of the most impor- 
tant challenges in contemporary cosmology. One usu- 
ally assumes the existence of a dark energy component 
which breaks the strong energy condition, and the sim- 
plest model that one can build is that of introducing a 
cosmological constant. Yet, in order to obtain its dom- 
ination today it has to be precisely set to an extremely 
small and so far unexplained value, one manifestation of 
the cosmological constant problem 0]. More generally, 
the fact that the dark energy overtakes dark matter at a 
recent epoch (z < 1) when we cosmologists are able to 
observe it is known as the coincidence problem. 

Allowing the dark energy to be dynamical may help to 
solve this issue. In this respect, a light scalar field, known 
as quintessence 0, 0, has been proposed. The field is 
meant to slow-roll down its potential, with its poten- 
tial energy acting analogously to that of early Universe 
inflation. A wide class of tracker models Q feature an 
attractor solution which roughly mimics the behaviour of 
the dominant component of the Universe, rendering the 
evolution of the field fairly independent of its initial con- 
ditions. Unfortunately, in order to obtain quintessence 
domination today, the parameters of the potentials so 
far discussed also need a fine-tuning, and so as yet those 
models have not led to a compelling resolution of the 
coincidence problem. 

More recently, models based on scalar fields with non- 
canonical kinetic energy 0], dubbed as k-essence 0,0) HI, 
have emerged. A subclass of models @ feature 
a tracker behaviour during radiation domination, and 
a cosmological-constant-like behaviour shortly after the 
transition to matter domination. As long as this transi- 
tion seems to occur generically for purely dynamical rea- 
sons, these models are claimed to solve the coincidence 
problem without fine-tuning. 

In this paper we will analyze the models given in the 
literature so far 00]. We will study the size of the basin 
of attraction of their tracker solutions and comment on 
the fine-tuning of the parameters. We will also look at the 
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behaviour of a general k-essence field close to singularities 
corresponding to a diverging sound speed. 

Throughout this article a prime denotes a derivative 
with respect to the argument of the function to which it 
is applied, and a dot denotes a derivative with respect to 
proper time. We assume 3/8nG = 1. 



II. K-ESSENCE FORMALISM 

In general k-essence is defined as a scalar field with 
non-canonical kinetic energy, but usually the models arc 
restricted to the Lagrangian 

£ k = K(cf>)p(X) , (1) 

where K(<f>) > and X = ±V ^V^. We note that 
the definition includes quintessence models (in this paper 
meaning scalar fields with canonical kinetic term) . Using 
the perfect fluid analogy, the pressure and the energy 
density are given by 

Pk (cf>,X) = K(0)p(X), (2) 

e k (0,A) = K(cf>)e(X), (3) 

where 

e(X) = 2Xp'(X)-p(X). (4) 

Following Refs. 0,0, we set K((/)) = l/4> 2 , define a new 
variable y = \j\fx and re-express p(X) as p(X{y)) = 
g(y)/y- In this case the pressure and the energy density 
become 



Pk(4>,y) = 

£k((f>,y) = 



g(y) 

cj) 2 y 



g'(y) 



(5) 
(6) 



We assume e^{4>,y) > 0, hence g'(y) < 0. The equation 
of state parameter and the effective sound speed are given 
by 
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We also assume w k (y) > — 1 and c 2 k (y) > which im- 
plies g" {y) > 0. As a result, g(y) must be a convex and 
decreasing function of y. 

From now on, we consider a flat Robertson- Walker 
Universe defined by the metric 



where 



ds 2 = -dt 2 +a 2 (t)dx 2 



(9) 



In this case, the Euler-Lagrange equation for the It- 
essence field is 



e'(X)4> + ZHp'(X)4> + §|| W = > 



(10) 



where H = a/ a. Then, if the Universe is filled with 
another fluid with energy density S{ and equation of state 
parameter wt constant one can find the following system 
of equations in terms of the independent variables y and 



yg"(y) 

3f2k (i - Ok) (wf - wk(y)) 
where N = ln(a/ao), f^k = £k/(£k + £f) and 



dy_ 
dN 
dQ k 
~d~N 



±] _Hg(y)-yg'(y)) . „ 



Here, we have assumed that <j> > 
and < fik < 1, the dynamics is co: 
trajectories in the y—Q^ plane. As 
system is well defined. As we can 
tionary point if w k (y s ) = W{ and r' 
in Ref. @ this is a stable point - 
sponds to a perfect tracking (i«k = 
fluid — if c 2 k (y s ) > w k (y s ). 



(11) 
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Therefore, for < y 
mpletely described by 
long as g"{y) ^ the 
see, (j/ s ,r2ks) is a sta- 
(y s ) — fiks- As shown 
- and therefore corre- 
= W{) of the dominant 



Ci 

C-2 



^-8g'(y c ) (r(y c ) - y^) 

Vc9"'{Vc) 
3{w k (y c )-l) 2(r(y c )-Vth) 



yc\/-Sg'(yc) 



(16) 
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If g'"{y c ) 7^ 0, Eq. I|15l) allows us to study the behaviour 
of the field close to the line y — y c , and unless fik = r 2 (y c ) 
the first term dominates. For the particular class of 
models we will analyze we have g'"(y c ) > 0, and there- 
fore if ilk < r2 (yc), hence C\ > 0, the solution moves 
away from the line, whereas if fl k > r 2 (y c ) the solu- 
tion ceases to exist (in that it hits the singularity y c ) 
within a finite time AN ~ 5y 2 /2\C\\ as it approaches 
the line (this has been checked numerically). We define 
S + = {(y c , ^k)|^k < r2 {yc)} as the segment of the singu- 
larity from which some trajectories spontaneously emerge 
and S_ = {(y c , fik)|^k > r2 (yc)} as the segment of the 
singularity on which some trajectories abruptly end. We 
can also use Eq. H15JI to determine the nature of the per- 
turbation Sy at the regular point So = (y c , ^kc). In that 
case the equation simplifies to give 



dSy 
dN 



:(wk(y c ) - l)5y. 



(18) 



III. SOUND SPEED DIVERGENCE 



Depending on the value of w(y c ) this solution either 
grows or decreases exponentially fast. In the cases we 
will study here, we have w k (y c ) < 1 and therefore by 
decreases. 

The existence of this singular behaviour means that 
a diverging sound speed leads to serious problems. In 
some situations it may well be possible to argue that 
the theory is valid up to a certain cut-off which excludes 
the singularity, but as we will show this is not the case 
for the cosmologically-realistic models proposed so far 0, 
8J. Instead, in these situations we must deal with the 
singular regions as we meet them. 



As we will see, for the particular class of models we will 
analyze, we have g"(y c ) = for some y c , which implies 
that i'(y c ) = and that the sound speed diverges at y c . 
In that case, from Eq. (jll|l we see that at y = y c there is 
a unique possible value Q k = Q kc given by the constraint 
equation 



^kc = r 2 (y c 



(14) 



This means that the phase space y-fik is cut into (at 
least) two parts separated by the line y = y c which is 
not allowed by the model except at So = (y c ,^kc)- To 
study the dynamics close to this line we use the expansion 
y = y c +Sy with Sy/y c <C 1, compute Eq. l(TT|) to 0(Sy 2 ) 
and find the equation 



dfry 
dN 
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C 2 5y-, 



(15) 



IV. ANALYSIS OF TWO MODELS 

K-essence models can possess many different attrac- 
tor solutions 0, Q, @ , especially trackers which perfectly 
mimic the dominant component of the Universe and at- 
tractors with negative equation of state leading to domi- 
nation of the field. By choosing an appropriate function 
g(y), it is possible to build a model with a certain number 
of attractor solutions which can feature some interest- 
ing dynamical properties. For instance, a stable tracker 
solution during radiation domination (R) could render 
the late-time evolution of the field fairly independent of 
its initial conditions. Then, the lack of such a solution 
during matter domination would force the field to reach 
another pseudo-attractor with w k ~ — 1 leading to the 
k-essence domination attractor (K) . This possibility was 
discussed in Refs. 0, @ and in that sense it is possible 
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TABLE I: This table shows the sign of some variables as a 
function of X for the two models analyzed in section llVl For 
each model the values of X c , X w and X* are given in the text. 
The symbol "oo" stands for a diverging value. 



to find a model which could solve the coincidence prob- 
lem: we cosmologists would observe the acceleration of 
the Universe today because we happen to appear soon 
after the onset of matter domination which leads to the 
formation of structures — and human beings. In this 
section we analyze the models proposed so far which are 
built in order to fulfill this principle @,@|. We study their 
sensitivity to initial conditions and comment on the fine- 
tuning involved. Throughout the discussion we will refer 
to and make use of the parameters first introduced in 
Refs. 0,0. 



A. Model 1 (Ref. [7]) 
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FIG. 1: Basin of attraction (shaded region) of the tracker so- 
lution R during radiation domination for the first model ana- 
lyzed in Section lTVl The points K and x denote the k-essence 
domination attractor and the saddle point respectively, de- 
fined in Ref. 8]. The solid lines are the limiting solutions 
which demarcate the basin of attraction and the long-dashed 
line separates the trajectories originating from and the ones 
originating from the singularity. The vertical dashed line is 
the singularity y — y c . So denotes the point (y c ,fikc)- 



Following the classification scheme introduced in 
Ref. this first model is of type (B r ). It is defined 
by 

p(X) = -2.01+2-/L + X+3x Kr 17 A 3 -Kr 24 A 4 . (19) 

As we will see, it appears to contain a number of prob- 
lematic issues which on the face of it contradict what is 
stated in Ref. 0. For example for some ranges of X 
we see that it is possible to have £k < 0, < — 1 and 
Cg k < 0. Moreover we find that c^ k and Wk diverge at 
X = X c w 1.6 x 10 7 and X = X w « 2.1 x 10 7 respec- 
tively. These properties are summarized in Tabled where 
we have also introduced X* sw 2.3 x 10 7 as the value of X 
for which the sound speed vanishes. From Eqs. J7J) and 
© , this also corresponds to the value where the equation 
of state parameter = — 1. In what follows we only an- 
alyze the dynamics for X < X c , which corresponds to 
y > y c rj 2.4 x 1CP 4 , since in this case > 0, % > — 1 
and Cg k > and also the solution which "solves" the co- 
incidence problem is in this region of the phase diagram. 

Since our concern over the dependency on initial con- 
ditions really only relates to the period of radiation dom- 
ination, we assume that the Universe is filled with radi- 
ation (wf = 1/3). We have run simulations in order to 
determine the size of the basin of attraction of the tracker 
R. In Fig. [l]we have plotted the phase diagram for the 
k-essence field during radiation domination. The solid 
and long-dashed lines are the limiting solutions which 
separate the different types of trajectories. The vertical 
dashed line is the singularity y — y c . As explained in 
section^] it can be divided in two parts: S+ (below So) 



and S_ (above So). All the trajectories contained within 
the long-dashed line originate from the same point and 
the rest originate from segment S+. The solid lines sepa- 
rate the trajectories ending on S_, the ones reaching the 
tracker solution R and the ones reaching the k-essence 
domination attractor K. (Having a look at the trajec- 
tories plotted in Fig. El can be useful for understanding 
the dynamics, although this figure describes the model 
we study next.) As we can see, the basin of attraction of 
R (shaded region) is very small and most of the solutions 
reach K. However as this is still during radiation domi- 
nation, they reach the k-essence domination attractor too 
early (typically after about 12 e-foldings, i.e. after an in- 
crease of the scale factor a(t) by about a factor 10 5 ) to 
be associated with the onset of matter-radiation equality 
and to be a candidate for dark energy. We also note that 
in a small region of the phase diagram the solutions cease 
to exist after a finite time as they reach S_, the part of 
singularity y — y c above So. 

Therefore, when analyzed in detail, this first model 
exhibits a rather different phase diagram structure from 
that sketched in Fig. 3 of Ref. to generically describe 
models of class (B r ). Unfortunately a complete com- 
parison is not possible as the precise equations used to 
generate that figure are not disclosed, so that we do not 
know to which k-essence model the figure corresponds 
(nor whether the model has an explicit lagrangian de- 
scription or is only of a more phenomenological nature). 
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FIG. 2: Basin of attraction (shaded region) of the tracker 
solution R during radiation domination for the second model 
analyzed in Section ITVI The explanations are the same as for 
Fig. Q S+ and S_ denote the parts of the singularity from 
where some trajectories originate and end respectively. A few 
trajectories have been plotted. 

B. Model 2 (Ref. [8]) 

This second model is defined by 

p(X) = -2.05 + 2^1 + /(X), (20) 

where 

f(X) = X - 1Q- & X 2 + 10~ 12 X 3 - IQ- 1& X A 

+ 10 -20 x 5_ 10 -24 x 6 /2 6_ ( 21 ) 

As in the previous example, this model is not always well 
denned. For X > X max ~ 6.3 x 10 5 we have l + f{X) < 
and therefore p(X) becomes ill-defined. Moreover, for 
X < X max , this model has the same problems as those 
of the model we have analyzed above. As before, its 
properties are summarized in Table [fl where this time 
X c » 3.9 x 10 5 , X w « 5.0 x 10 5 and X* m 5.3 x 10 5 . 

As in the first example, this is a model of type (B r ) 
and for y > y c sa 1.6 x 10~ 3 its phase space has a simi- 
lar structure to that shown in Fig. ^ In Fig. [2 we have 
plotted the basin of attraction of the tracker solution R 
during radiation domination for this second model. In 
order to help in understanding the dynamics a few tra- 
jectories have been added. Again, we see that for most of 
the initial conditions the field does not reach the tracker, 
but instead it reaches either an early k-essence domina- 
tion solution K or the singularity on S_. Therefore, this 
second model suffers the same fine-tuning as that of the 
first example. 

V. DISCUSSION 

As introduced in Refs. II], the idea of k-essence 
seems appealing. Current models of quintessence suf- 



fer in general because of the fine-tuning of the potential 
parameters to account for the fact that the field has only 
recently started dominating. K-essence was introduced 
as an extension of quintessence models by taking into ac- 
count non-canonical kinetic terms. In a subclass of these 
models a tracker behaviour occurs during radiation 
domination, and a cosmological-constant-like behaviour 
shortly after the transition to matter domination. Since 
the k-essence field seems to change its behaviour gener- 
ically for purely dynamical reasons, these models could 
be claimed to solve the coincidence problem without fine- 
tuning of the initial conditions. In that sense they are 
more natural than quintessence models which rely on a 
different, non-obvious, scale for the transition to occur. 

In this short paper, we have addressed in a bit more de- 
tail the question over the nature of this attractor solution 
during radiation domination which is meant to avoid a 
fine-tuning of the initial conditions. We have numerically 
solved the evolution equations for the k-essence fields and 
written the results in terms of the physically motivated 
parameters y and £lk, following Refs. @|- The key 
result we have found is that the basin of attraction for 
the tracker solution appears to be very small compared 
to the basin of attraction for the k-essence domination 
solution and therefore it cannot be seen as equivalent to 
that of quintessence models. In other words, for almost 
all initial conditions, the system would evolve rapidly into 
k-essence domination. It would have done so way before 
matter-radiation equality and in that sense, the required 
behaviour of k-essence can only hold for a specific subset 
of initial conditions. If these turn out to be an impor- 
tant set of conditions then k-essence can be thought of 
as providing an elegant way of obtaining the acceleration 
we see today. However, if there is no particular reason 
for choosing such initial conditions, then we believe that 
k-essence suffers from the same fine-tuning issues that 
plague quintessence models. 

With regard to this issue, in Refs. 0,0 the authors ar- 
gue that the basin of attraction of the radiation attractor 
of the two models studied in this paper are compatible 
with cquipartition. This particular initial condition may 
be the case should the k-essence or quintessence field be 
associated with one of the many fields produced at the 
end of a period of inflation. However, it need not be the 
case, and if it was not, equipartition would not help in 
choosing initial conditions. Indeed, equipartition is gen- 
erally applied to systems where particle production oc- 
curs after the decay of the inflaton, whereas in the case of 
vacuum energy the initial condition may well have been 
set by early Universe physics (see for example Ref. 0). In 
any case, both basins of attraction are so small that their 
main parts do not overlap and therefore there remains a 
fine-tuning issue even when assuming equipartition. 

We have also shown that the models proposed so far 
feature a singularity associated with a diverging sound 
speed. This leads to problems like the sudden disappear- 
ance of some field trajectories. The presence of this sin- 
gularity could be avoided by adding an extra term to the 
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functions given by Eqs. and JUJ, but our attempts 
to do so have led to the creation of a second radiation 
attractor which goes against the goal of having a model 
for which the late-time behaviour is independent of the 
initial conditions. Another way of solving this problem 
would be to argue that the theory is valid up to a cer- 
tain cut-off which excludes the singularity, but as long as 
the radiation tracker is very close to the singularity this 
argument cannot be applied. 

In this paper we have dealt with the specific la- 
grangians given in Refs. 0, H| at face value. Those la- 
grangians have a complicated form which is not appeal- 
ing, and it is disappointing that so far it has not proved 
possible to produce much simpler models, as one would 
expect to be able to if the desired triggering behaviour re- 
ally is generic to k-essence. The particular nature of the 
models makes it hard to assess how difficult it is to build 
a working example, for instance quite how much fine- 
tuning is required to establish the time needed between 



matter radiation equality and k-essence domination, and 
how easily one can avoid the presence of many different 
attractors during radiation domination. Unfortunately, 
we are not aware of any particular particle physics moti- 
vated models which would deliver the lagrangians of the 
two models, and, as with quintessence, we believe that 
a realistic model of k-essence remains a challenge which 
has to be met. 
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